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Resume 

Apres un bref apergu permettant de situer notre travail, nous pro- 
posons une nouvelle voie pour aborder la programmation lineaire en 
proposant un algorithme elabore a partir d'une idee simple qui per- 
met d'obtenir une solution aussi approchee que voulu par translation 
dichotomique d'un hyperplan de l'espace des solutions. 

Mots cles : Algorithme, Programmation Lineaire, polyedre, polytope, hyper- 
plan, dichotomie, optimisation, approximation. 

Abstract 

After a short course in order to situate our work, we propose a new 
way to study linear programming and we give a proposal of algorithm 
to solve linear programming from a basic idea which allow to obtain 
an approached solution with desired accuracy. For this we use some 
dichotomic translations of an hyperplan in the solutions hyperspace. 
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1 Introduction 

Le probleme de progr animation lineaire s'est rapidement impose des qu'on 
a voulu planifier un tant soit peu les activites economiques ou autres. C'est 
ainsi que des les annees 1939, les necessites de la planification sovietique 
conduisent Kantorovitch et Tolstoi a proposer une solution au probleme, (voir 
[Kan39[ rTol39j ) inspiree plus ou moins des travaux de Joseph Fourier (1768- 
1830). En 1951 Dantzig publie les resultats de Kantorovitch- Tolstoi sous 
forme d'algorithme executable sur ordinateur (voir |Dan51j ) et lui donne le 
nom de simplex. Cet algorithme fait le tour du monde sans aucun concur- 
rent jusqu'en 1979 avec l'apparition de l'algorithme dit de I'ellipsoide (voir 
[Kha79] ) . qui est en temps polynomial, ce qui, dans le pire des cas (degeneres- 
cence), n'est pas le cas du simplex. Parait ensuite un autre algorithme polyno- 
mial base sur une autre idee, celui de Karmarkar (voir [Kar84aj et [Kar84bj ). 
puis en 1989 l'algorithme de Murty & Chang (voir [CM89j ). Ce dernier algo- 
rithme etant ameliore en 2006 (voir [Mur06bt IMur06a] ) . Ces algorithmes, sim- 
plex, ellipsoi'de, Karmarkar, Murty sont souvent caracterises comme etant, de 
frontiere pour le simplex, exterieur pour l'algorithme de I'ellipsoide, interieur 
pour les algorithmes de Karmarkar et de Murty. D'autres aussi proposent 
une approche un peu differente (voir par exemple [NPL091 IPL09j ). Dans le 
present article, nous avons voulu presenter une methode, qu'on peut qualifier, 
d' exterieure-interieure, ou encore alternante basee sur une idee geometrique 
de dichotomie de l'espace. II s'agit de considerer (dans le cas d'une maximisa- 
tion) un hyperplan particulier de la famille de ceux dermis par la fonctionnelle 
a optimiser, tantot exterieur au polytope des contraintes, tantot le recoupant, 
et de se rapprocher par translation controlee, autant qu'on le souhaite de la 
solution optimale par un argument de dichotomie. On obtient ainsi une solu- 
tion approchee aussi precise que voulu, ce qui est essentiel dans le domaine 
pratique, 1' optimum sur le mo dele n' etant quasiment jamais 1' exact optimum 
dans la realite modelisee. 
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2 Position du probleme 

Dans toute la suite E designe un R— espace-vectoriel ; R etant l'ensemble 
des nombres reels. 

On appelle probleme de programmation lineaire, tout probleme pouvant 
s'ecrire sous la forme : 



max Ym=i °i x i Q G R; n G N* 
s/c ( ^ i=1 aijXj ^ b * a, b G R; j = 1 : m, i = 1 : n (1) 

On rappelle qu'un polyedre convexe de est un ensemble P de la forme 

P = {x G lR n : Ax < 6 }, (2) 

ou A : E — > R est une application lineaire (m G N; si m = 0, P = P,) 
b G R m et Finegalite Ax < b se lit composante par composante dans R m 
(Ax)i < bi pour tout i G {1, • • • ,m}. Si l'ensemble se presente avec des 
egalites lineaires, Cx =, d on pourra se ramener sous la forme ([2]) en les 
remplgant par deux inegalites opposees Cx < d et — Cx < —d. Un polytope 
est un polyedre convexe et borne. 

Geometriquement, un polyedre est done l'intersection d'un nombre fini 
de demi-espaces de P. Si P est de dimension finie, il n'y a pas de restriction 
a supposer que P = W 1 et que A est une matrice m x n (il suffit de donner 
une base de P). Dans certaines circonstances, par exemple en optimisation 
lineaire, il est plus avantageux de representer un polyedre de R n sous la forme 
dite standard suivante : 



P = {x G R n Ax = b, x > }. (3) 

II n'y a aucune perte de generalite dans cette representation. Tout polyedre 
de la forme ([2]) se represente sous la forme ([3]) en introduisant des variables 
d'ecart. 

Remarque 2.1 Les representations et |3|) sont dites duales, car elles 
font intervenir des applications lineaires (elements du dual de E). 

On considere le polytope D (voir figure) suppose contenir l'origine O, sinon, 
on peut toujours introduire les variables d'ecart de sorte que l'origine O G D. 
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Separation de I'espace en deux parties par un hyperplan dont I'une contient 
le polytope. 

3 Quelques rappels en optimisation 

Soient / une forme lineaire sur E, et a G R, A, B deux sous ensembles de 

E. 

n 

f(X) = kjXj, n e N*; kj e R 

3=1 

Definition 3.1 ?7n sous-ensemble C de E est dit convexe si Vx, y G C et 
A G [0, 1] on a Xx + (1 — A)y G C , On note 

[x,y] = {Xx + (l-\)y,0<X<l} 

Definition 3.2 On appelle hyperplan H d 'equation [f = a], V ensemble definit 
par 

H={X G R N , tel que f(X) = a.} (4) 
H est ferme <^==^> / est continue. 
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3.1 Separation des ensembles convexes 

Un outil essentiel en analyse convexe est le theoreme de Hann-Banach sur 
la separation des ensembles convexes. On supposera que l'espace vectoriel E 
est de dimension finie. On peut toujours le munir d'un produit scalaire note 

<•,•>■ 

La separation de deux ensembles convexes se fait geometriquement dans 
E en utilisant un hyperplan affine H, e'est-a-dire un ensemble de la forme 

H = {x G E : = a }; 

ou £ G E est non nul et a G E. On dit que cet hyperplan separe deux 
convexes C\ et Ci si l'on a 

Vxi G Ci, Vi 2 G C 2 ■ (£,xi) < a < {£,x 2 ) 

Ceci sera certainement le cas s'il existe un £ non nul dans E tel que 

sup (£,xi) < inf (£,x 2 ) 

On dit que cet hyperplan separe strictement ces deux convexes s'il existe 
deux scalaires a>i et a 2 tels que ctx < a < a 2 et 

Va?i G Ci, Vi 2 6 C 2 : (£, x x ) < a x < a 2 < (£, x 2 ) 

Ceci sera certainement le cas s'il existe un £ (necessairement non nul) dans 
E tel que 

sup < inf (f,x 2 ) 

La figure ci dessous illustre cette notion. Ici, on a utilise le produit eucli- 
dien sur IR 2 . 
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Separation de deux convexes 



3.2 Theoremes de Hann-Banach 

Nous allons donner ici quelques resultats classiques en optimisation convexe 
et dont les demonstrations peuvent etre trouvees dans [Gil06j . et aussi dans 
|Ach84] pour ce qui concerne les polyedres. Ces resultats sont connus dans la 
litterature sous le nom de formes geometriques du theoreme de Hann-Banach 
ou de separation des ensembles convexes. Le premier affirme que Ton peut 
separer strictement deux ensembles convexes disjoints, si Fun est ferme et 
l'autre est compact. Le second exprime que Ton peut separer (non stricte- 
ment cette fois) deux ensembles convexes quelconques (en dimension finie). 

Theoreme 3.1 (Hann-Banach I) SoitE un espace vectoriel de dimension 
finie, muni d'un produit scalaire note (.,.). Soient C\ et Ci deux convexes 
disjoints non vides de E, tels que R Cf° = {0}. Alors, on peut separer 
C\ et C*2 strictement : il existe un vecteur £ G E tel que 

sup < inf (£,x 2 ) 

Corollaire 3.1 Soient C\ et C2 deux convexes non vides disjoints d'un es- 
pace vectoriel E muni d'un produit scalaie note (., .), on suppose I 'un ferme 
et l'autre compact, alors on peut separer C\ et C2 strictement. 
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Theoreme 3.2 (Hann-Banach II) Soit E un espace vectoriel de dimen- 
sion finie, muni d'un produit scalaire note (.,.). Soient C\ et C2 deux convexes 
disjoints non vides de E. Alors, il existe un vecteur £ £ E non nul tel que 



Definition 3.3 Soient C C E un convexe et H un hyperplan ferine tel que 
3u £ C fl H, on peut supposer que H = {/ = a} f £ £" etant le dual de 
E, et a £ R. C < {/ > a} cm C C {/ < a} on dit que H supporte A en 
u ou que u est un point de support de C . 

Corollaire 3.2 1- Soit C un convexe ferme non vide Vw £ dC (bord de C ), 
u est un point de support de C ; Si C est un ensemble convexe d'un espace 
vectoriel topologique, par tout point frontiere passe un hyperplan d'appui. 
2- Soit C un ensemble convexe ferme, Alors 



p est V ensemble des demi-espaces fermes contenant C . 

4 Algorithme propose 

Dans ce paragraphe, nous proposons un algorithme d'optimisation per- 
mettant de resoudre le probleme de programmation lineaire donne sous forme 



Soit H l'hyperplan defini par (j3J) et K le polytope defini par l'ensemble 
des contraintes de (pQ). En utilisant les theoremes de separation, on peut 
toujours separer Fespace en deux parties par l'hyperplan de sorte qu'une 
partie contienne le polytope K des contraintes. 

Trois possibilites apparaissent alors : 

- soit l'hyperplan H ne coupe pas le polytope K des contraintes et dans le 
cas d'une maximisation comme dans le PL de la formule[T]et alors H est 
<k au-dela > des solutions realisables, e'est-a-dire qu'il n'en peut conte- 
nir aucune et il convient alors de choisir un autre hyperplan H / / H 
situe entre l'origine du repere (i.e. le point 0) et H ; 

- soit l'hyperplan H coupe le polytope des contraintes et il faut alors in- 
clure H dans les contraintes et choisir un autre hyperplan de separation 
de l'espace H*//H situe dans le demi-espace ne contenant pas 0; 



sup (£,xi) < inf (£,x 2 ) 
xieCi X2GC2 
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- soit le polytope des contraintes est vide. 
Si H n K = 0, alors, en reprenant la notation de la definition 13.21 on definit 
un nouvel hyperplan H' — {/ = ^} par exemple; il va exister alors deux 
valeurs de a, a SQ et a Sl , avec s_i et s G N* telles que l'hyperplan H SQ = 
{/ = a So } fl if 7^ alors que H s ± — {/ = a S l } fl if = 0. On pose alors : 

H s+1 = {f = s+i = } 

Puis si i^s+i Pi ^ = on P ose s -i = s +i i 
Sinon, si H s+1 f] K ^ 0, on pose so = s+i. 
Dans cet algorithme, l'hyperplan choisi est defini par 

n 

H = {x e R n I ^CiXi = a} (5) 
i=i 

oix Y^i=i °i x i represente la fonction objectif du probleme de programmation 
lineaire. 

Remarque 4.1 A chaque pas de V algorithme, le nouvel hyperplan se deduit 
du precedent par translation. En effet, on cherche les vecteurs qui engendrent 
Vespaee H defini par |3]) puis on cherche un vecteur unitaire 17" qui represente 
le pas vectoriel de la translation 



4.1 Un sous-probleme cle 

Dans cet algorithme, une difficulte reside dans Identification de la va- 
cuite -ou la non-vacuite- de l'intersection de l'hyperplan avec le polyedre des 
contraintes : 

#p|/i=0? (6) 

Au plan theorique, Particle de Murty |Mur06b] resoud le probleme polyno- 
mialement. Du point de vue pratique, dans les cas de non-degenerescence, le 
probleme peut etre resolu par l'utilisation de la methode dite " du grand M" . 



5 Conclusion et Perspectives 

Cette etude ouvre une voie originale en matiere de programmation 
lineaire qui n'est inspiree d'aucune des methodes existantes mais qui utilise 

1. Cette etude nous a ete facilitee par l'amitie et l'aide bibliographique de Pham Canh 
Duong du CNST a Hano'i. 
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toutefois la methode dit du point interieur de Murty pour la resolution d'un 
sous probleme cle. II reste a examiner en quoi cette methode est utilisable en 
programmation lineaire en nombres entiers et en programmation lineaire pa- 
rametrique. La methode peut aussi s'etendre a la programmation non lineaire 
dans quelques cas particuliers. 
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